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ABSTRACT 

We investigate the ultraviolet-to-optical spectral energy distributions (SEDs) of 17 active 
galactic nuclei (AGNs) using quasi-simultaneous spectrophotometry spanning 900-9000A (rest 
frame). We employ data from the Far Ultraviolet Spectroscopic Explorer (FUSE), the Hubble 
Space Telescope (HST), and the 2.1-meter telescope at Kitt Peak National Observatory (KPNO). 
Taking advantage of the short-wavelength coverage, we are able to study the so-called "big blue 
bump," the region where the energy output peaks, in detail. Most objects exhibit a spectral 
break around llOOA. Although this result is formally associated with large uncertainty for some 
objects, there is strong evidence in the data that the far-ultraviolet spectral region is below the 
extrapolation of the near- ultraviolet-optical slope, indicating a spectral break around llOOA. We 
compare the behavior of our sample to those of non-LTE thin-disk models covering a range in 
black-hole mass, Eddington ratio, disk inclination, and other parameters. The distribution of 
ultraviolet-optical spectral indices redward of the break, and far-ultraviolet indices shortward of 
the break, are in rough agreement with the models. However, we do not see a correlation between 
the far-ultraviolet spectral index and the black hole mass, as seen in some accretion disk models. 
We argue that the observed spectral break is intrinsic to AGNs, although intrinsic reddening as 
well as Comptonization can strongly affect the far-ultraviolet spectral index. We make our data 
available online in digital format. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: general — ultraviolet: general 
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1. INTRODUCTION 

The spectral energy distribution (SED) of AGNs contains a significant feature in the ultraviolet (UV) 
to optical region, known as "the big blue bump." This feature is thought to be thermal emission from an 
optically thick accretion disk feeding a massive black hole (e.g., Shields 1978; Malkan & Sargent 1982), and 
it has been argued that its energy peak lies in the unobserved extreme ultraviolet (EUV, ~100-912A) region 
(e.g., Mathews & Ferland 1987). To determine the shape and the peak of the big blue bump is of particular 
importance, since this provides critical information on the structure and condition of the inner-most region 
in AGNs as well as on the ionizing flux that powers the emission lines. 

Zheng et al. (1997) constructed composite AGN spectra from HST spectra, and found that the UV- 
optical power-law continuum breaks at around 1000A. This was later confirmed in a similar HST composite 
using a larger sample (Telfer et al. 2002). Laor et al. (1997), also using composite spectra, found that the 
soft X-ray continuum matches up with the extrapolation of the HST composite, consistent with the idea 
that the UV bump actually peaks in far-ultraviolet (FUV, 912-2000A), rather than the EUV. However, such 
a spectral break is not seen in the recent FUSE composite spectra for low-redshift AGNs (Scott et al. 2004), 
although some individual sources do show a spectral break. Scott et al. (2004) noted that this FUSE sample 
has a median luminosity one order of magnitude lower than that of the HST sample (Zheng et al. 1997), 
and it is possible the break wavelength is luminosity dependent. It is important to realize that composite 
spectra represent a complex average of many individual spectra. Examination of the individual spectra in 
a composite may lead to greater insights, and they can be used to understand the overall characteristics of 
a composite in a more physical way. For instance, the spectrum of 3C273 shows a break near the Lyman 
limit (Kriss et al. 1999), and reasonable fits to the spectrum with accretion disk models have been reported 
(Kriss et al. 1999; Blaes et al. 2001). 

It is well accepted that AGNs are powered by accretion onto massive black holes. Many accretion 
disk models have been built to predict the AGN continuum and compare with observed spectra (e.g., Sun 
& Malkan 1989; Laor & Netzer 1989; Laor 1990; Czerny & Zbyszewska 1991; Hubeny et al. 2000, 2001). 
Since the AGN local environment is not well known, and the models are still relatively very simple, it is 
not possible for the models to predict the detailed features in the observed spectra. However, large-scale 
features in the AGN SEDs can be reproduced by disk models. Simple models produce a large continuum 
discontinuity at the Lyman limit, but it is not often seen in real objects (e.g., Koratkar, Kinney, & Bohlin 
1992). Such a feature could be smeared out by relativistic smearing and Comptonization (e.g, Hsu & Blaes 
1998; Kriss et al. 1999; Blaes et al. 2001; Hubeny et al. 2001). This results in a bump or a spectral break, 
instead of an edge, in the vicinity of Lyman limit, which resembles the spectral breaks seen in some observed 
AGN spectra. Hubeny et al. (2000) constructed geometrically thin accretion disk models with non-LTE 
atmospheres, including a full treatment of general relativistic effects in the disk structure. We choose this 
specific model and compare the spectral properties in the FUV-optical region for a sample of 17 low-redshift 
AGNs with the model expectations. 

We characterize the spectral continuum with broken power-laws. Unless noted, the power-law indices 
we use through this paper are all a\ (f\ cx A" A ) except that the soft X-ray spectral index a x is a u (f u cx 
v au ). It is easy to convert between a\ and <x v (a\ + ot v = —2). For cosmology we choose A = 0, H = 
75 km s _1 Mpc _1 , and q a = 0.5. Since the objects in this sample are at low rcdshifts, the results in the paper 
are not affected by the choice of cosmology. 
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2. SAMPLE AND DATA 

The FUSE AGN program (Kriss 2000) has surveyed more than 100 of the UV-brightest AGNs, of which 
about 20 were also observed in an HST spectral snapshot survey during 1999-2000. The FUSE observations 
were scheduled as close in time as possible with the HST snapshot observations. Ground-based optical 
spectra were also obtained during the same period at KPNO. We excluded a few objects because of the 
lack of an optical spectrum (NGC 3783, low declination), strong host galaxy contamination (NGC 3516), or 
strong variability (NGC 5548, also no simultaneous HST spectrum) . As a result, we have compiled a sample 
of 17 AGNs, with quasi-simultaneous spectrophotometry covering rest wavelength from 900-9000A. This is 
a heterogeneous sample with low redshift (z < 0.5). 

Table 1 lists the basic parameters and observation log for the sample. For more than half of the 
objects, the FUSE, HST, and KPNO spectra were obtained within a few months; all but 5 objects (3C351, 
IRAS F07546+3928, Mrk509, PG0947+396, and PG1 100+770) were observed in the three bands within a 
year. In the case of PG1100+770, we failed to observe the red part of the spectrum at KPNO in Feb 2000, 
and have instead used archival (1993) spectrophotometry obtained with the 2.7 meter telescope at McDonald 
Observatory that was consistent with the 2000 epoch blue KPNO spectrum. For 3C351 and Mrk509, archival 
HST data are used, but the FUSE and optical observations were essentially simultaneous, and can be used 
to constrain the flux level of the archival HST spectra when necessary. Our HST and optical observations 
of IRAS F07546+3928 were close in time, but the contemporaneous FUSE observation had now signal due 
to wrong pointing of the telescope. We therefore use a FUSE observation from the FUSE archive obtained 
later in time for this object. 



2.1. Optical Spectra 

All the optical spectra were obtained with the 2.1m telescope at KPNO except for a few noted in Table 1. 
A wide slit of 6" was used to ensure that all the light from the AGNs is included. Two spectra were obtained 
for each object, covering the observed wavelengths from ^3180-6000A and ^5600-9000A, with resolution of 
~9A and <~12A, respectively. 

We used standard packages in IRAF to reduce the optical data. We paid special attention in subtracting 
the host galaxy contribution when extracting one-dimensional spectra. A low-order polynomial was fit across 
the dispersion direction within the extracting aperture to represent the host galaxy contribution and sky 
background. The galaxy contribution was clear in many of the lowest luminosity objects, and appeared to 
be well subtracted from the final AGN spectra. We estimate that the host galaxy contamination is less than 
5% in all cases. 

Wavelength calibration was done by using comparison spectra, and absolute flux calibration was achieved 
by using the standard star spectra obtained on the same night when it was photometric. At least one of 
the two observations for each object was done under photometric conditions. The two spectra were then 
combined in the observed frame. Usually, the flux calibration of the two spectra agree very well; in a few 
cases when they do not match in the overlap region, we scale one spectrum to match the one with better 
flux calibration. 



Table 1. Sample and Observation Log 



Observation Date Dataset ID CalFUSE 



Object 


Other Name 




z a 


E(B-V) b 


FUSE 




HST 




Opt. 


Bluc 


Opt. 


Red 


FUSE ID 


HST ID 


Ven 




3C273 


PG1226+023 





1576 





.021 


2000- 


■ 04- 


■23 


2000- 


■03- 


■16 


2000- 


■ 02- 


■25 


2000- 


02-26 


P1013501 


O5G045.IZQ , 


O5G045K0Q 


v2 


.0 


5 


3C351 


PG1704+608 





.3730 





023 


1999- 


■ 10- 


■17 


1992- 


■ 02- 


■15 


1999- 


■ 10- 


■09 


1990- 


09-20 


Q1060101 


Y0VM0103T 


Y0RV0G04T C 


v2 


2 


.2 


4C + 34.47 


B2 17214-34 





2055 





037 


2000- 


■06- 


■09 


2000- 


■06- 


■25 


1999- 


■ 10- 


■09 


2000- 


02-27 


P1073501 


O5G077Y2Q . 


O5G077Y3Q 


v2 


.2 


2 


IRAS F07546+3928 


MS 0754.6 + 3928 





0953 





066 


2002- 


■ 02- 


■ 11 


2000- 


■ 01- 


■29 


1999- 


■ 10- 


■09 


1999- 


10-08 


S6011801 


O5G018LPQ 


, O5G018LQQ 


v2 


.1 


7 


MRK290 


PG1534+580 





0303 





015 


2000- 


■03- 


■16 


2000- 


■06- 


■02 


2000- 


■ 02- 


■25 


2000- 


02-27 


P1072901 


O5G070TLQ 


, O5G070TMQ 


vl 


.8 


7 


MRK304 


PG2214+139 





0657 





073 


2000- 


■07- 


■16 


2000- 


■06- 


■19 


1999- 


■ 10- 


■07 


1999- 


10-08 


P1073901 


O5G082ANQ 


, O5G082AOQ 


v2 


.0 




MRK506 







0428 





.031 


2000- 


■06- 


■08 


2000- 


■06- 


■ 24 


1999- 


■ 10- 


■09 


2000- 


02-27 


P1073401 


O5G076UEQ 


, O5G076UFQ 


v2 


.0 




MRK509 







.0345 





057 


1999- 


■ 11- 


■06 


1992- 


■06- 


■ 21 


1999- 


■ 12- 


■ 11 


1999- 


-12-11 


X0170102 d 


Y0YA0302T , 


Y0YA0305T e 


v2 


1 


.7 


NGC3516 







0883 





042 


2000- 


■ 04- 


■17 


2000- 


■ 04- 


■20 












P1110404 


O5G032T7Q , 


O5G032T8Q 


v2 


.1 


7 


NGC3783 







0097 





.119 


2000- 


■ 02- 


■02 


2000- 


■05- 


■17 












P1013301 


O5G039LBQ 


O5G039LCQ 


vl 


.8 


7 


PG0052+251 







1544 





.047 


1999- 


■ 10- 


■03 


1999- 


■ 10- 


■01 


1999- 


■ 10- 


■07 


1999- 


10-08 


P1070101 


O5G003NAQ 


, O5G003NBQ 


v2 





.5 


PG0947+396 







2057 





.019 


2001- 


■ 01- 


■06 


2000- 


■06- 


■15 


1999- 


■ 10- 


■09 


2000- 


-02-27 


A0600101 


O5G023NPQ 


, O5G023NQQ 


v2 


2 


.2 


PG0953+414 







2338 





.013 


1999- 


■ 12- 


■30 


2000- 


■ 02- 


■05 


1999- 


■ 10- 


■09 


2000- 


02-26 


P1012202 


O5G024NBQ 


, O5G024NCQ 


v2 


2 


2 


PG1100 + 772 


3C249.1 





3114 





.034 


2000- 


■ 01- 


-20 


2000- 


■ 01- 


■31 


1993-0 


5 f 


1993-05 f 


P1071601 


O5G030D4Q 


O5G030D5Q 


v2 


.2 


.2 


PG1259 + 593 







.4769 





008 


2000- 


■ 02- 


■25 


2000- 


■ 02- 


■09 


2000- 


■ 02- 


■25 


2000- 


-02-26 


P1080101 


O5G047IJQ , 


O5G047IKQ 


v2 


.2 


.2 


PG1322+659 







1684 





019 


2000- 


■05- 


■08 


2000- 


■06- 


■18 


2000- 


■ 02- 


25 


2000- 


-02-26 


A0600808 


O5G052WXQ 


, O5G052WYQ 


v2 


.2 


.2 


PG1351+640 







0882 





.020 


2000- 


■ 01- 


■18 


1999- 


■ lo- 


■28 


2000- 


■ 02- 


■25 


2000- 


-02-26 


P1072501 


O5G054KQQ 


, O5G054KRQ 


v2 


.1 


.7 


PG2349-014 


PKS 2349-10 





1740 





.027 


2000- 


■06- 


■25 


1999- 


os- 


■27 


1999- 


■ 10- 


■07 


1999- 


10-08 


P1074201 


O5G088N3Q . 


O5G088N4Q 


v2 


.2 


2 


TON951 


PG0844+349 





0643 





037 


2000- 


■ 02- 


-20 


1999- 


■ 10- 


21 


2000- 


02- 


28 


2000- 


02-27 


P1012002 


O5G020QBQ 


, O5G020QCQ 


vl 


.8 


7 



a Measurcd from the optical data in this study (§4.4). 
b From NED based on Schlcgel, Finkbeincr. & Davis (1998). 
c Obscrved on 1991-10-22, Y0RV0G03T is also used. 
d X0170101 is also used. 

C Y0YA0303T and Y0YA0304T arc also used. 

^Obtained with 2.7m telescope at McDonald Observatory. 
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2.2. Near Ultraviolet Spectra 

Near-UV (NUV) spectra covering a wavelength range of ~1150-3180A were obtained from our HST 
spectroscopic snapshot survey. Observations were made with STIS in slitless mode, to minimize target 
acquisition overheads. A guide star acquisition assured good imaging and hence full spectral resolution 
(~lA) on these point sources. To save more overhead time, we skipped the standard wavelength calibration 
observations that are used to determine the zero-point correction for the wavelength scale. Instead, we used 
the Galactic absorption lines clearly visible in each spectrum to determine the proper zero-point correction. 
To apply this to the spectra, we use an iterative process. We run the CALSTIS pipeline with no zero-point 
correction, and then measure the wavelengths of the Galactic lines. We then calculate the required correction, 
assuming that the Galactic lines are at their rest wavelengths, and manually edit the header information in 
the CALSTIS pipeline at the point that the wavelength zero-point correction is made. We then measure the 
wavelengths again in the extracted spectra. Usually one iteration is sufficient, and more than two are never 
required. The resulting errors are typically less than 1 A. 

The spectra taken with G140L and G230L were then combined. When the two spectra do not match in 
the wavelength overlap region, we scaled the G230L spectrum by a uniform factor, which typically resulted 
in a correction less than a few percent. 

2.3. Far Ultraviolet Spectra 

FUSE spectra covering observed wavelengths of 905-1 187A have a high resolution of ^0.05 A. Newer 
versions of the standard FUSE calibration pipeline CalFUSE were used to process the raw FUSE/FW data 
(see Sahnow et al. 2000). Spectra were extracted after background subtraction with updated background 
models and subtraction algorithms. The wavelength and flux were then calibrated. The version of CalFUSE 
used for different objects are also listed in Table 1 for reference. 

Among all the instrumental effects in the FUSE data, the most prominent is the "worm," a dark stripe 
of decreased flux in the spectra running in the dispersion direction. The flux loss can be as much as 50% in 
the longer wavelengths, where we need to connect with HST spectra. To correct for the effects of the "worm" , 
we used data from the two independent LiF channels covering the 1100-1180 A wavelength range. The worm 
is most often present in channel LiFlb. After confirming this by inspection, we form a ratio of the spectra 
obtained in the two independent channels. Assuming that the data in channel LiF2a are uncorrupted, we 
fit a low-order spline to the ratio to derive a correction curve for LiFlb, and apply this to the LiFlb data. 
This removes the "worm" while preserving the high-spectral resolution and statistical independence of the 
LiFlb data. More information can be found at the FUSE webpage 1 . 

2.4. Soft X-ray data 

We do not have simultaneous X-ray observations. Instead, we collected available soft X-ray spectral 
indices and fluxes (0.1-2.4 keV, Table 3) from the literature (Brinkmann, Yuan, & Siebert 1997; Pfefferkorn, 
Boiler, & Rafanelli 2001). These data were obtained from both R OS A T All-Sky Survey (RASS) and public 
PSPC observations. Assuming Galactic absorption, the power-law X-ray photon indices were estimated from 



1 http://fuse.pha.jhu.edu/analysis/calfusc.html 
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the two hardness ratios given by the Standard Analysis Software System (Brinkmann, Yuan, & Siebert 1997) 
or from a spectral fit (Pfcfferkorn, Boiler, & Rafanelli 2001). X-ray fluxes between 0.1-2.4keV were also 
obtained. We used the photon index and the X-ray flux to calculate the flux density at 1 keV and the X-ray 
spectral index a x , and used them in the correlation analysis (§4.5). 

We did not find X-ray information for PG 1259+593, and the X-ray flux of Mrk 304 is too low to derive 
a spectral index (Brinkmann, Yuan, & Siebert 1997). There are Chandra observations of NGC 3783 (e.g., 
Kaspi ct al. 2002), but the data show strong variability and warm absorbers that complicate a consistent 
comparison with the spectral indices derived from the ROSAT observations for other objects. 

3. CONSTRUCTION OF SPECTRAL ENERGY DISTRIBUTIONS 

The very strong geocoronal emission lines, i.e. airglow emission (e.g., OiA1304, LyaA1216, Ly/3A1026 
etc.), were first removed by hand from both the FUSE and HST spectra. We made no attempt to remove 
numerous narrow Galactic molecular hydrogen absorption lines. The FUSE, HST, and optical spectra were 
then combined. The overlap regions between the spectra arc usually 20-100A, large enough for us to assess 
the agreement of flux level in the spectra due to flux calibration, source variability, and "worm" correction for 
FUSE spectra. The overlap region between the HST and optical spectra is relatively small, and sometimes 
there is a gap of 5-40A. However, we found that this does not prevent us from matching the continuum 
levels on both sides of the gap, because a flux change of 5% at the gap would be very obvious by checking 
the overall continuum trend on both sides of the gap. 

The flux densities of the optical and HST spectra usually agree very well within 3%, and we do not scale 
the spectra to match in this case because the uncertainty in the scaling factor is at a similar level. When 
the flux density of the FUSE spectrum does not agree with the above, we scale it to match using a constant 
scaling factor determined from the overlap regions between FUSE and HST spectra (Table 2). In the cases 
when the FUSE, HST, and optical spectra do not agree with each other, we scale the FUSE and optical 
spectra to match the HST spectrum, because based on our knowledge, the flux calibration of the HST STIS 
spectra is usually more reliable. Most flux differences are slight (< 20%), and corrections in the 20-40% 
range affect less than a quarter of the sample. Since the spectral shape, rather than the absolute flux, is the 
most important in constructing SEDs, this scaling process should not significantly affect our study of SEDs 
if problems in the flux calibration are the cause (e.g., because of clouds). 

However, based on our knowledge of the flux calibrations, we believe the change of the continuum levels 
in different wavebands is most likely due to source variability. Although we tried to obtain quasi-simultaneous 
spectra for each object, some observations were actually separated by a few months, and the continuum level 
could vary by a factor of 20% or more, and it is empirically true that most AGNs show spectral shape 
changes with flux changes. We list the observing time gaps and scaling factors applied for the spectra in 
Table 2 as a guide to assess this problem. 

We corrected for Galactic reddening with an empirical mean extinction law (Cardelli, Clayton, & Mathis 
1989, CCM), assuming Ry = Ay / 'E(B-V) = 3.1, a typical value for the diffuse interstellar medium. 
E{B — V) is obtained from NED 2 based on the dust map created by Schlegcl, Finkbcincr, & Davis (1998). 
Since the Cardelli et al. (1989) extinction curve 's UV cutoff is at lOOOA, we extrapolate the curve down to 



2 NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and Space Administration, http://nedwww.ipac.caltech.edu/ 
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Tablc 2. Observing Time Gaps and Scaling Factors 



Object 


Observing Time Gap (days) a 




Scaling Factor 




FUSE 


HST 


opt. blue 


opt. red 


FUSE 


HST 


opt. bluc D 


opt. rcd D 


3C273 


38 





-20 


-19 


1.03 


1 


1 


0.94 


3C351 





-2801 


-8 


-2288 


1.10 


1 


0.76 


scaled 


4C+34.47 


-16 





-260 


-119 


1.13 


1 


1 


1 


IR07546+3928 


744 





-112 


-113 


1.12 


1 


1 


0.95 


MRK290 


-78 





-98 


-96 


1.33 


1 


1 


0.73 


MRK304 


27 





-256 


-255 


0.99° 


1 


1 


0.95 


MRK506 


-16 





-259 


-118 


1.09 


1 


l d 


l d 


MRK509 





-2694 


35 


35 


1.01° 


1 


1 


1.03 c 


PG0052+251 


2 





6 


7 


1.01° 


1 


1 


1 


PG0947+396 


205 





-250 


-109 


0.85 


1 


1 


scaled 


PG0953+414 


-37 





-119 


21 


0.67 


1 


1 


scaled 


PG1100+772 


-11 





~-2450 


~-2450 


0.91 


1 


1.10 


1.10 


PG1259+593 


16 





16 


17 


0.77 


1 


1 


1 


PG1322+659 


-41 





-114 


-113 


1.04 


1 


0.79 


1 


PG1351+640 


82 





120 


121 


1.00 


1 


1.67 


1.50 


PG2349-014 


303 





41 


42 


0.44 


1 


1 


1 


TON951 


122 





130 


129 


0.82 


1 


1 


1 



a Rclative to HST observing time except for 3C351 and Mrk509, for which HST archival data are used and the 
time gap is relative to FUSE observing time. 

b Whcn the optical blue spectra match the HST spectra within a few percent, the same level as the uncertainty, 
we do not apply a scaling (scaling factor— 1). Optical red spectra are often not photometric. Scaling the red to 
match the blue is not a problem since they were usually obtained within 1-2 days. 

c No scaling is applied when combining the spectra since the scaling factor is small. 

d Data were photometric, but host galaxy contribution is large. We removed the galaxy contribution until 
everything matches properly. The correction of galaxy is good to 10%. 
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900A for our FUV spectra by using the same formula. This extrapolation of CCM law to shorter wavelength 
has been shown to be compatible with recent FUSE data (e.g., Hutchings & Giasson 2001). 

We finally applied the rcdshift correction and brought both the wavelength (in vacuum) and flux density 
(fx) of the spectra to rest frame. The redshifts were determined using measurements of [Oni] A5007 in our 
optical data (§4.4). 

The final SEDs of this sample are shown in Figure 1. All the spectra of the 19 AGNs (including 
NGC 3516 and NGC 3783, but not NGC 5548) are available at http://physics.uwyo.edu/agn/. These 
include the unsealed individual FUSE, HST, and optical spectra as well as the combined SEDs. The special 
treatment of wavelength calibration in the HST spectra is not a problem for using the data for other studies, 
e.g., emission- line analyses, — the wavelength calibrations are only slightly less accurate than standard STIS 
observations (0.5-1.0 pixels rms vs. 0.1-0.2 pixels), and this has no measurable effect on the flux calibration. 

4. SED ANALYSIS 

4.1. SED Measurements 

For many objects, the FUV and NUV-optical spectral regions have different slopes (Fig. 1). After careful 
examination, we decided to use three power-laws to fit the entire FUV-optical continuum. 

A power-law with spectral index of chfuv is fitted to the FUV spectral region ( <1100A) on a case-by- 
case basis. We first exclude the obvious emission-line regions (Ly7A973, CiiiA977, NiiiA991, Ly/3A1026, 
OviA1034, SiivA1062, SiivA1073, and HenA1084, etc.) Strong ISM Lyman series absorption lines and 
possible AGN O vi and Lya absorption features are also excluded. We make sure the excluded regions are 
wide enough by visual inspection so that the line wings are also excluded. We then fit the remaining data 
points with a power-law, and iteratively reject points beyond ±2a and the points next to the rejected points. 
This process removes the numerous narrow Galactic H 2 absorption lines and possible residual emission 
features from the fitting. Finally, we visually inspect the fitting results and make sure that the power-law 
goes through the apparent continuum regions in the FUV spectra. We repeat this process to estimate 
upper and lower limits for apuv by intentionally including some emission or absorption features until the 
fitted power-law obviously deviates from the spectra (visual inspection). We therefore obtain conservative 
uncertainties for the best-fit power-law. The fitted cxfuv an d its uncertainties are listed in Table 3. If there 
are possible residual weak blended Galactic absorption lines or emission features, they are only comparable 
with the noise level, and the uncertainty in the fitted spectra indices caused by these weak features is much 
smaller than the above estimated uncertainties. 

A single power-law cannot fit the entire NUV-optical region in many objects. This was noticed before, 
for example, in the Sloan Digital Sky Survey (SDSS) composite spectra (Vanden Berk et al. 2001). Clean 
continuum regions are also hard to find in the AGN NUV-optical spectra due to the large number of broad 
emission lines and blends, including the "small blue bump" from ^2000-4000A — the blend of Fell emission 
and Balmer continuum. We therefore have to define, for this sample, some common, narrow continuum 
windows where there seem to be no emission lines: 1144-1157A, 1348-1358A, 4200-4230A, 5600-5648A, 
6198-6215A, and 6820-6920A. 

The NUV-optical spectra index auvo m ^1200-5500A, and red optical spectra index aored in ^5500- 
9000A are each obtained by fitting a power-law to a pair of selected continuum windows, requiring that all 
emission features are above the fitted power-laws in the corresponding regions. The continuum windows 
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Rest Wavelength (A) 



Fig. 1. — FUV to optical SEDs and fitted power-laws for different regions. Uncertainties in the FUV power- 
law fitting are also shown. The vertical dotted lines indicate the continuum windows used for fitting auvo 
(-1200-5500 A) and a 0re d (-5500-9000A) (Table 3). FUSE spectra have been rebinned to a resolution of 
0.5 A for display purpose. 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 2. — Distribution of the break wavelength, indicating a spectral break near HOOA for most objects. 
The largest X break is from IRAS F07546+3928 (§4.2) and PG2349-014. Mrk506 does not show a spectral 
break and is not included. 



used for each object are listed in Table 3 and marked in Figure 1. Since the continuum windows are very 
narrow, this fitting process is more like defining each power-law with two points. The power laws cannot 
be treated as the true continua of the spectra, but the spectral indices provide information on the overall 
continuum slopes. 

The UV bump can now be characterized with two power-laws (a broken power-law) with spectral indices 
of Q.FUV and auvo, and a break wavelength Xbreak, which is defined by the intersection of the two power-laws. 



4.2. Spectral Break in SEDs 

As can be seen in Figure 1, an extrapolation of the NUV-optical power-law does not match the FUV 
continuum in most objects. For the possible exceptions, Mrk 290, Mrk 506, Mrk 509, PG0947+396, and Ton 
951, the extrapolated NUV-optical power-law falls within the bounds of the errors for our fits to the FUV 
continuum. Thus, for 12 out of 17 objects, we see a break in the spectral index to a steeper value when 
comparing the NUV-optical to the FUV continuum. 

The break wavelength is calculated as the intersection point of the two power-laws. The distribution of 
the break wavelength is shown in Figure 2, where Xbreak peaks near 1100A and spans 800-1600A. However, 
since the calculated break wavelength is very sensitive to small changes in a fuv or auvo when the difference 
between cxfuv and auvo is small, it has a large uncertainty for some objects (Table 3). 

We have also compared our SEDs with the soft X-ray spectral indices in Figure 3 in a similar way as 
in Laor et al. (1997), except that we also have FUV data. For more than half our objects, the soft X-ray 
spectral indices appear to match up reasonably with the extrapolation of the FUV continuum. This directly 
confirms the finding by Laor et al. (1997) and Zheng et al. (1997) in composite spectra that the peak of the 
big blue bump lies in the FUV region. 



Three objects, IRAS F07546+3928, NGC 3516, and PG1351+640, show strongly suppressed NUV-FUV 
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continua. These objects also show intrinsic absorption features, possibly suggesting the existence of dust 
associated with the absorbers. (See Zheng et al. (2001) for the case of PG 1351+640.) We will discuss the 
reddening effect more in §6. 

4.3. Bolometric Luminosity 

Bolometric luminosity (Lb i) is one of the most fundamental parameters for understanding the black 
hole accretion in quasars, however, it has not been easy to obtain for quasars in general, because they emit 
significant power over a large part of the electromagnetic spectrum. 

Elvis ct al. (1994) built SEDs for a sample of 47 quasars and were able to obtain the bolometric lumi- 
nosity by integrating over the SEDs from radio to X-ray wavelengths. They also determined the bolometric 
correction factors for a few monochromatic luminosities, e.g., Lb i ~ 13ALa(5400A), based on their SEDs. 
Different empirical correction factors have been determined and used to estimate Lb i m previous studies 
(e.g., Sanders et al. 1989; Laor & Draine 1993; Wandel, Peterson, & Malkan 1999). Recently, many studies 
use the prescription of Kaspi et al. (2000) to estimate quasar black hole masses and bolometric luminosity 
L Bo i = 9AL A (5100A). 

Since we have broad spectral coverage from the FUV to optical, we are able to obtain an accurate 
luminosity for this region by integrating over the power-laws we measure. To estimate the bolometric 
luminosity, we need to include X-ray and IR regions. We extend our FUV power-law continuum to 700A, 
and then use a power-law to connect with the soft X-ray luminosity at 0.2 keV. 

The case for the infrared region is more complicated. There is an IR bump around 10^m in quasar 
SEDs, which also contributes a significant amount of energy. After averaging the mean SEDs for radio-loud 
and radio-quiet objects from Elvis et al. (1994), we fit two power-laws to characterize this bump and obtain 
a\ = —0.69 for 1-lO^m and a\ = —1.65 for 10-100/zm. These two power-laws form a peak at 11.45/im, 
roughly corresponding to the peak in the mean SEDs. We scale this power-law IR bump to match the 
extrapolation of the fitted NUV-optical continuum at 1/im for each object. 

We estimate the bolometric luminosity between 2keV and 100/zm by integrating this set of power-laws 
over this region. We have also obtained the luminosities for individual wavebands. Table 4 lists the results. 

Since we do not use actual measurements in the IR (few exist for our sample objects), we scale the 
IR bump in two ways so that we have a range of the estimated IR luminosity. In case A, we match the 
IR bump with the extrapolation of the red optical power-law of aored at l^m; in case B, we ignore aored 
and use NUV-optical power-law of auvo for the entire optical region, and match the IR bump with the 
extrapolation of auvo at 1/im. Therefore, we have two integral luminosities for FUV-optical Lpuvo (700A- 
1/zm), two estimated luminosities of the IR bump Lir (l-100^m), and two estimates of Lb i (2 keV-100yum). 
These two cases give consistent results for Lpuvo, but result in a big difference, up to a factor of 2.5 in 
Lin, when auvo an d aored are significantly different. In any case, Lju is comparable with Lpuvo and 
contributes significantly to L Boi ■ Table 5 shows the comparison between Lir and L fuvo ■ The distribution 
of Lin/ Lpjjvo has a large dispersion, in general agreement with the data of Elvis et al. (1994). 

Figure 4 compares our integral Lb i with the bolometric luminosity estimated using the empirical 
formula from a monochromatic optical luminosity, Lb i = 9ALa(5100A). The ratios of Lb o ;/9ALa(5100A) 
are listed in Table 4. (The ALa(5100A) is measured in a local continuum. See §4.4). There exists a strong 
correlation between the integral Lb i and 9ALa(5100A), but our results are 30-70% (<~ 0.1 — 0.2 dex) larger 
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Table 3. Spectral Indexes 





FUV 


Cont. Windows' 3 


1200-5500A 


5500- 


-9000A 






0.1-2.4 kcV 


Object 


/l000 a a FUV 


a b c d c f 


/iooo a a uvo 


/iooo a 


a Ored 




x d 


/lkcv C a x° 



3C273 


53 


.70 


-0 


3C351 


3 


.44 


-0 


4C + 34.47 


1 


.95 





IRAS F07546... 


1 


.95 


-0 


MRK290 


1 


.99 


-0 


MRK304 


5 


.78 





MRK506 


2 


.15 


-1 


MRK509 


14 


50 


-1 


PG0052 + 251 


4 


.41 


-0 


PG0947+396 


1 


.92 


-1 


PG0953+414 


7 


.43 


-1 


PG1100+772 


4 


.38 


-0 


PG1259+593 


4 


.17 


-0 


PG1322 + 659 


2 


.26 


-1 


PG1351 + 640 


1 


.96 





PG2349-014 


1 


.57 


-0 


TON951 


5 


.74 


-1 



0.64 
0.28 
0.00 
0.22 



+ 0.24 
0.11 
+ 0.30 
0.28 
+ 0.21 
0.61 
+ 0.35 
0.54 



-0.66 
+ 0.52 
-0.01 
+ 0.16 
'"-0.38 
,, + 0.91 
-0.23 
+ 0.45 
0.04 
+ 0.37 
0.02 
+ 0.13 
-0.04 
+ 0.44 
-0.14 
+ 0.11 
0.12 
+ 0.31 
0.14 
+ 0.51 
0.20 
+ 0.31 
0.22 
+ 0.33 
0.27 



0.81 
1.33 



0.64 
1.06 
0.60 
0.70 
1.05 



60.00 
4.69 
2.39 
2.50 
2.50 
8.52 
1.92 

15.30 
4.84 
1.84 
8.35 
4.99 
5.72 
2.51 
2.35 
2.32 
5.86 



-1.74 
-1.35 
-1.39 
-0.73 
-1.55 
-1.44 
-1.11 
-1.54 
-1.67 
-1.52 
-1.81 
-1.73 
-1.80 
-1.66 
-0.57 
-1.51 
-1.37 



0.97 
5.99 
0.48 
7.15 
1.29 
4.32 
2.19 
0.60 
2.91 
0.57 

0.62 
7.80 
0.49 
5.12 



-0.84 
-1.22 
-0.44 
-1.30 
-0.88 
-0.80 
-1.17 
-0.81 
-1.19 
-0.50 

-0.85 
-1.28 
-0.60 
-1.24 



1.10 
1.07 
1.39 
0.51 
0.94 
2.18 
0.02 
0.31 
0.86 
0.19 
0.76 
1.26 
1.16 
0.60 
1.17 
0.81 
0.32 



+ 
-69 
T +57 



1105 
1027^ 

ie2o+i«» 

1281+", 
,+49 
-43 



1155 T 



1193^ 



1 9m +358 
i2ui -167 



+ 78 
-48 
+ 323 
-87 
, , rr + 39 

+ 10 
-72 
+ 1 

-93 

1191 + 36 
11J -196 
+ 101 



1111 
813 



1114 
1313 



1173 
1616 
1066 



-62 
+ 292 
284 
+ 880 



283.0 
3.1 

66.3 
6.9 

30.1 

26.6 
386.0 
46.2 
14.9 
17.0 
12.5 

17.8 
4.2 

31.2 
3.2 



-1.11 
-1.31 
-1.29 
-2.16 
-1.32 



+ 0.01 
0.01 
+ 0.05 
0.05 
+ 0.06 
0.06 
+ 0.38 
0.38 
+ 0.13 
0.13 



-1 20+ - 09 

-1 61 + ' 03 
° -0.03 

-1 49+ ' 02 
-0.04 
„ + 0.18 
-0.15 
, + 0.03 
-0.05 
_, «. R +0.08 
-^"O-O.IO 



-1.18 
-1.43^ 



75 + 0.04 
-0.03 
-1 43+0-06 
1 '* i, -0.06 
_1 44+0-12 
1 '**-0.12 
4+0.11 
-0.12 



-1.75^ 



-1.54^ 



a ^1000 — fitted rest frame continuum flux density at 1000A for corresponding regions (10 14 erg s 1 cm 2 A 1 ). 

b Continuum windows used for fitting NUV-optical region. a:1144-1157; b:1348-1358; c:4200-4230; d:5600-5648; c:6198-6215: f:6820-6920. For Mrk506, 1444-1458A 
is used instead of a; for PG2349-014, 5550-5650A and 7200-7300A arc used instead of c and f. 

c Ao = a FUV - a uvo 

^Break wavelength for ctpjjy and &JJVO- Th° errors arc calculated solely from the errors of ctpjjy since the uncertainty of ct^j V Q is negligible based on the way 
it is measured (§4.1). 

6 ^lkev — f\ at 1 kcV (10~ 14 erg s~ 1 cm -2 A -1 ). Calculated from integrated flux between 0.1-2.4keV and ct x (Ju « u Clx ) from Brinkmann, Yuan, & Siebert 
(1997) and Pfefferkorn, Boiler, & Rafanclli (2001, for Mrk290 and Mrk506 only), assuming power-law. 



Table 4. Integral Bolometric and Individual Waveband Luminosities 



Object 


L_ 


K 


Lxfuv 




Lfuvo 




L 


IR 












I 














case 


A a 


case 


B a 


case A a 


case 


ga 


case 


A a 


Ra 


case 


ga 


Rf3 ° 


3C273 


45. 


.51 


46.22 


46 


.51 


46 


.51 


46.42 


46. 


41 


46 


.94 


1 


.48 


46 


.94 


1.47 


3C351 


44. 


23 


45.44 






46 


.04 




46. 


20 










46 


.50 


1.03 


4C+34.47 


45. 


13 


45.30 


45 


.42 


45. 


.41 


45.72 


45. 


56 


46 


.07 


1 


.62 


46 


.00 


1.39 


IR07546+3928 


43. 


.83 


44.54 


45 


.09 


45. 


.11 


45.53 


45. 


64 


45 


.72 





.98 


45 


.80 


1.17 


MRK290 


43. 


23 


43.64 


43 


.86 


43 


.81 


44.26 


43. 


86 


44 


.52 


1 


.83 


44 


.34 


1.20 


MRK304 








45 


.02 


45 


.01 


45.22 


45. 


.16 


45 


.46 


1 


.08 


45 


.42 


0.99 


MRK506 


43. 


44 


43.95 


44 


.20 


44 


.19 


44.54 


44 


48 


44 


.81 


1 


.39 


44 


.78 


1.29 


MRK509 


44. 


.53 


44.87 


44 


.75 


44. 


.74 


44.96 


44. 


77 


45 


.44 


2 


.28 


45 


.38 


1.99 


PG0052+251 


44. 


.81 


45.33 


45 


.43 


45. 


42 


45.53 


45. 


37 


45 


.98 


1 


.83 


45 


.93 


1.62 


PG0947+396 


44. 


45 


45.11 


45 


.31 


45 


.29 


45.54 


45. 


32 


45 


.86 


1 


.71 


45 


.76 


1.36 


PG0953+414 


44. 


67 


45.66 


45 


.94 


45 


.94 


45.94 


45. 


78 


46 


.41 


1 


.72 


46 


.35 


1.53 


PG1100+772 


44. 


78 


45.63 


45 


.95 


45 


93 


46.13 


45 


85 


46 


.47 


1 


.79 


46 


.36 


1.37 


PG1259+593 C 












46 


.24 




46 


13 










46 


.57 


1.09 


PG1322+659 


44. 


.54 


45.16 


45 


.22 


45. 


.21 


45.36 


45 


16 


45 


.80 


1 


.93 


45 


.73 


1.66 


PG1351+640 


43. 


30 


44.13 


45 


.08 


45. 


11 


45.52 


45 


71 


45 


.68 





.83 


45 


.83 


1.16 


PG2349-014 


44. 


.71 


45.07 


45 


.22 


45. 


.20 


45.54 


45. 


30 


45 


.87 


1 


.82 


45 


.77 


1.45 


TON951 


42. 


95 


44.39 


44 


.92 


44 


90 


45.12 


45. 


05 


45 


.41 


1 


.26 


45 


.37 


1.15 



Note. — Values are the logarithm of luminosity in units of ergs s 1 . Lx- 2— 0.2keV; Lxfuv'- 0.2kcV-700A; Lfuvo'- 
700A-l^m; L IR : l^m-lOO^m: L Bo i- 2kcV-100/xm. 

a Case A: IR bump is scaled to match aor-ed at 1/xm; Case B: IR bump is scaled to match auvo at l/i m - a Or*ed is ignored 
and auvo ls used for the entire optical region. 

|3 R j 4,Rb: ratios of the integral LboI^-o 9ALa(5100A) for case A and B, respectively. 

c Luoi docs not include Lx and Lxfuv- 
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- PG0052+251 
- PG0947+396 
- PG0953+414 
- PG1 100+772 
- PG1259+593 
PG1322+659 
- PG1351+640 
- PC2349-014 
T0N951 




Fig. 3. — FUV-optical SEDs and the soft X-ray spectral indices. The dotted lines are drawn to connect FUV 
spectra and X-ray spectral slopes, and the dashed lines indicate the extrapolation of fitted FUV power-law. 
More than half of the objects have an X-ray spectral slope roughly matching the extrapolation of the FUV 
continuum slope. Also note the large UV suppression of IRAS F07546+3928, NGC 3516, and PG1351+640. 
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in general and seems to agree better with the correction factor determined by Elvis et al. (1994). These 
imply that the bolometric correction factor may be larger than 9, and more like 13 for using AL A (5100A). 

Although we obtain accurate FUV-to-optical luminosities from our SEDs, a few sources contribute to 
the uncertainties of the integral bolometric luminosity. First, Elvis et al. (1994) warned about the use of 
mean SEDs, because the shapes of SEDs in individual objects have large dispersion. Second, our scaling of 
the IR bump may cause large uncertainties since we use simple extrapolations of the NUV-optical and red 
optical power-laws. We have seen a difference of 30% (0.1 dcx) in Lb i for case A and B. Third, we have 
no data for Lxfuv, an d the cutoff wavelength of 700A for Lfc/v is chosen subjectively. Finally, the X-ray 
data were not obtained simultaneously with our FUV-optical spectra, and the X-ray variability is another 
source of uncertainty. These should be kept in mind when interpreting our integral LboI- 

For simplicity and easy comparison with other work, we will use -Lb o /=9ALa(5100A) for the bolometric 
luminosity in the rest of the paper. This choice does not change any statistical significance in analyses 
involving Lb i, and it is easy to compare with our integral bolometric luminosity by using Table 4 to correct 
for individual objects. 



4.4. Estimation of Black Hole Mass 

We have estimated the black hole mass and accretion rate using a recently developed method based on 
reverberation mapping of the broad-line region (BLR) and on the assumption of virial motion (Kaspi et al. 
2000), 

M BH = Rblr v 2 /G. (1) 

H/3 is used for estimating the velocity dispersion, v = V3/2 FWHM(H/3), and Rblr is the size of the broad 
line region and can be estimated empirically from reverberation mapping studies (Kaspi et al. 2000), 

Rblr = 32.91?;' 

We use the bolometric luminosity Lb a i — 9ALa(5100A). Given the black hole mass and bolometric luminosity, 
we can also estimate the Eddington ratio, L/Lsdd- 

We measure the FWHM(H/3) by fitting the H/3 region with the IRAF task specfit (Kriss 1994). A local 
power-law continuum, the [Oni] lines, and HenA4686 are also fitted together with H/3 (Fig. 5). We use 
a broad and a narrow Gaussian component to fit the H/3 broad line, and allow a relative wavelength shift 



Table 5. Statistics of Lib,/Lfuvo 





Median 


Mean 


Min 


Max 


Case A 


1.61 


1.76 ±0.63 


0.78 


2.78 


Case B 


1.20 


1.49 ±0.91 


0.70 


4.00 


Elvis et al. (1994) a 


1.17 


1.56 ±1.23 


0.58 


5.89 



'^Lfuvo is between 0.1-1/im. Data are obtained from their 
Table 15 for 34 objects, for which L(l-10^m) and L(10-100Mm) 
are not given as upper limits. 



AL A (5100A) 
10 44 erg s- 1 



light days. 



(2) 
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between the two components to account for the H/3 asymmetry. A narrow-line-region (NLR) component of 
H/3 is also introduced, but it is often negligible. Both the width and the wavelength of this NLR component 
are tied with those of [Oin] A5007. The intensity ratio of [Oiii] A5007/A4959 is assumed to be 3:1 based on 
their statistical weights. A single Gaussian profile is fitted for each of [Oiii] AA4959,5007. 

Broad Fe n emission blends are often strong in this region, especially when [O III] is weak. In order to 
remove the Fen contamination, we use an Fen template from Boroson & Green (1992). The strength of the 
Fe II template is free to vary in the fitting process, and it is also broadened to be consistent with FWHM(H/3) 
for each object. 

We calculate the final FWHM(H/3) from the fitted H/3 model profiles. Specifically, we exclude the NLR 
component, and use the other two fitted Gaussian components, taking into account the relative wavelength 
shift between these two components. The uncertainty of FWHM(H/3) is less than 10%. 

The rest-frame reference is defined with [O in] A5006.8, and the fitted wavelength of [Oiii] A5007 is used 
to calculate the redshift. The uncertainty of the rcdshift is < 0.0002 for all objects except for PG1259+592, 
which has very weak [Oin], and has a redshift uncertainty of 0.002. 

Table 6 lists the fitting results, together with calculated Mbh and L/LEdd etc. Six objects in our sample 
also have black hole masses derived from reverberation mapping studies (Kaspi et al 2000). They are also 
listed in the table. Our results are usually larger, because we have excluded the H/3 NLR component when 
obtaining the FWHM(H/3). Therefore we have a larger FWHM(H/3) and hence a larger Mbh- As was also 
noticed by Boroson (2002) and Vestergaard (2002), 3C351 (i.e., PG 1704+608) is an extreme case, in which 
the very narrow H/3 component on top of the broad H/3 emission is identified as an NLR component with 
the same width as [Oiii] (690km s _1 ). It is excluded from calculating FWHM(H/3) in our study, but not in 
Kaspi et al. (2000) 3 . This results in a huge difference in FWHM(H/3) and hence in M BH . 

Another uncertainty of the FWHM(H/3) comes from the uncertainty of the fitted local continuum level. 
Assuming a single Gaussian profile, if the continuum is lowered by 10% of the line peak, the estimated 
FWHM will increase by ~ 7%, and the calculated black hole mass from FWHM(H/3) and the continuum 
luminosity based on Eq.l and 2 will increase by 7%. 

On the other hand, sometimes we had to scale the optical spectra to match the HST NUV spectra 
due to likely source variability. A scaling of 30% can change the estimated black hole mass by 20%. The 
AGN intrinsic variability can also change the emission-line profile, and hence the estimated black hole mass. 
Different approaches for spectral measurements in different studies can lead to large discrepancies in reported 
masses. Due to all the above reasons, the estimated black hole mass from different studies can easily differ 
by a factor of a few, in our case, a maximum factor of 5 for 3C273 (excluding the extreme case of 3C351). 
In fact, Vestergaard (2002) compared the black hole masses estimated from reverberation mapping and 
from single-epoch optical observation, and concluded that they agree within factors of 3, 6, and 10 with 
probabilities of 80%, 90%, and 95%, respectively. However, if one keeps consistency in measurements and 
calculation for a sample, the relative uncertainty in the estimated black hole masses within the sample should 
be much smaller. 



3 Their measured mean FWHM(H/3) is 890 km s 1 , and it is 400 km s 1 from rms spectrum. The rms spectra in the 
reverberation mapping studies are not completely free of constant narrow components (Peterson et al. 1998). 
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Fig. 4. — Comparison between integral Lb i and 9ALa(5100A). The solid line indicates Lb i — 1-5 x 
9ALa(5100A), it is not a least square fit. 







Table 6. 


Black Hole Mass and Eddington 


Ratio 






Object 


z 


FWHM(H/3) 


/ A (5100A) a 


AL A (5100A) b 


M BH 


L/L E dd 


M B h (rcv) c 






(km s" 1 ) 


(org s _1 cm~ 2 A" 1 ) 


log(org s _1 ) 


(1O 8 M0) 




(1O 8 M0) 


3C273 


0.1576 


4115 


3.41E-14 


45.82 


11.50 


0.41 


2.35 


3C351 


0.3730 


9760 


4.11E-15 


45.54 


41.18 


0.06 


0.075 


4C+34.47 


0.2055 


3520 


2.64E-15 


44.91 


1.94 


0.30 




IR07546+3928 


0.0953 


2965 


7.89E-15 


44.78 


1.12 


0.39 




MRK290 


0.0303 


5505 


2.38E-15 


43.31 


0.36 


0.04 




MRK304 


0.0657 


5570 


7.96E-15 


44.48 


2.43 


0.09 




MRK506 


0.0428 


5520 


3.15E-15 


43.72 


0.70 


0.05 




MRK509 


0.0345 


3630 


1.22E-14 


44.13 


0.59 


0.16 


0.92 


PG0052+251 


0.1544 


5465 


3.12E-15 


44.77 


3.73 


0.11 


3.02 


PG0947+396 


0.2057 


3810 


1.53E-15 


44.68 


1.57 


0.22 




PG0953+414 


0.2338 


3155 


4.23E-15 


45.22 


2.57 


0.46 


1.64 


PG1100+772 


0.3114 


9300 


2.94E-15 


45.27 


24.20 


0.06 




PG1259+593 


0.4769 


3615 


3.11E-15 


45.58 


6.03 


0.45 




PG1322+659 


0.1684 


3030 


1.67E-15 


44.56 


0.82 


0.32 




PG1351+640 


0.0882 


2840 


9.69E-15 


44.81 


1.08 


0.43 


0.30 


PG2349-014 


0.1740 


5900 


1.98E-15 


44.66 


3.64 


0.09 




TON951 


0.0643 


2390 


6.19E-15 


44.36 


0.37 


0.45 





a Rest- frame flux density at 5100A. 

b Assuming zero cosmological constant, Hq — 75 km s _1 Mpc _1 , and qo — 0.5, same as the cosmology used in Kaspi et al. 
(2000). 

c Black hole mass measured from reverberation mapping studies (Kaspi ct al. 2000). 
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Fig. 5. — Model fitting to the H/3 region. The dotted lines show the data, and the thick solid lines the fitting 
results. Also shown are the local continua, Fen template, and individual components of H/3, [Oin], and 
He II A4686. 
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4.5. Correlation Analyses 

It is natural to think that the properties of the UV bump, thought to arise from the accretion disk, arc 
governed by the AGN fundamental parameters, such as Mbh and L/ LEdd- For example, standard thin disk 
models predict low disk temperatures for high Mbh and/or low L/LEdd (Shakura & Sunyaev 1973) and 
therefore longer break wavelengths. We looked for such correlations in our sample, but we do not find any 
significant correlation between the UV bump properties we have measured (ocfuv, ctuvo, Aa, and Xbreak) 
and Mbh, or L/LEdd- Table 7 lists the Pearson correlation coefficients for selected parameters. A principal 
component analysis has also been performed on these parameters, and no hidden correlations are revealed. 
We also group the objects into subsamples based on L/Lsdd (or Mbh), but there is still no evidence of 
a correlation within subsamples (see Fig. 6 for an example). However, Scott et al. (2004) found a weak 
correlation at 96% confidence level between apuv and Mbh m their subsample of 21 AGNs. There are only 
4 objects in common between these two studies. The inconsistency between the results most likely arises 
from the small size of the samples. In any case, the fact that there are no or weak correlations suggests that 
Mbh and L/Lsdd are not the only parameters underlying the observed properties of the UV bump. Other 
factors, such as the disk inclination, intrinsic reddening or other unidentified parameters must also play an 
important role (§5). We also note that a large spread over Mbh and L/Lsdd within the small samples may 
wash out any correlation with spectral index, but our sample is too small to address this. 

A strong correlation (r = 0.82, p = 0.0001) seen in Table 7 is between Aa (= apuv — otuvo) and apuv- 
This is simply because the distribution of the NUV-optical spectra index auvo is relatively narrow, while 
otFUV spans a wide range. Figure 7 shows this clearly. Most objects have auvo between —2 and —1 with a 
median value of —1.52, but the distribution of czfuv is broader (Fig. 7b). 

There seems to be an anti-correlation between a x and Xbreak, but there are only 14 objects with values 
of both a x and Xbreak {r = — 0.67, p = 0.01). Since Xbreak has very large uncertainty, this correlation should 
not be treated seriously. 

We also see a correlation (r = 0.64, p = 0.01, for 15 objects) between E(B-V) and a x , but this is largely 
due to an outlier, IRAS F07546+3928, with the softest a x = —2.16. Without this outlier, the correlation 
disappears (p — 0.23). After careful checking, we find no other correlations that are created or destroyed by 
outliers. 



Table 7. Pearson Correlation Coefficients (r) 





\og(M BH ) 


L/L Edd 


a FU v 


auvo 


Act 


Xbreak 


E(B-V) a x 


log(M Bff ) 


1.00 














L/L Edd 


-0.16 


1.00 












a F uv 


0.23 


-0.06 


1.00 










auvo 


-0.33 


0.12 


0.46 


1.00 








Aa 


0.47 


-0.14 


0.82 


-0.13 


1.00 






^break 


0.02 


0.06 


0.34 


0.39 


0.14 


1.00 




E(B-V) 


-0.16 


-0.24 


0.31 


0.30 


0.15 


0.28 


1.00 




0.23 


-0.24 


-0.09 


-0.31 


0.15 


-0.67 


-0.64 1.00 



Note. — 17 objects are used except for a x (15 objects) and Xbreak (16 objects). The chance 
probability (p) of 1%, 2%, and 5% corresponds to a correlation coefficient of 0.61, 0.56, and 0.48, re- 
spectively, for 17 objects. Large correlation coefficients marked in bold face do not reveal significant 
physical correlations (see §4.5 for detail). 
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Log M BH (M ) 



Fig. 6. — No correlation between o>fuv and Mbh for the total sample or for subsamples grouped by L/LEdd- 
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Fig. 7. — Distributions of auvo and olfuv of the sample (solid line), all models (dotted line, §5) and 
the selected models (dashed line, 1O 8 M < Mbh < 4 x 1O 9 M , 0.03 < L/L E dd < 0.3). Note that no 
selected models are as red as auvo > — 1- Two reddest objects (auvo > — 1) are IRAS F07546+3928 and 
PG 1351+640, both of which show evidence of intrinsic dust reddening. 
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5. COMPARISON WITH THIN-DISK MODELS 

As mentioned above, the expected correlations between the UV bump and Mbh or L/LEdd may be 
mitigated by the small sample size, a large variation in these parameters within the sample, and by other 
parameters affecting the spectrum of the accretion disk. We used the thin disk model developed by Hubeny 
ct al. (2000) to investigate this. 

The models are constructed for a non-LTE disk with 5 free parameters: black hole mass, mass accretion 
rate (M), viscosity parameter a V i SCO , black hole spin, and inclination angle cosi. The total spectrum of a 
disk is integrated over the individual annuli, taking into account the inclination angle. We have chosen a 
grid of models with a maximally rotating Kerr black hole with possible values of Mbh between 0.125 x 10 9 
and 32 x 10 9 M Q , M between 2^ 12 M Q yr~ 1 and 64M Q yr~ 1 , a visco = 0.01 or 0.1, and cosi between 0.01 and 
0.99. The maximum Eddington ratio is limited to L/LEdd ~ 0.3. Above this value, the model disk becomes 
geometrically thick and thus no longer self-consistent. 

In order to make statistical comparisons between theory and observation, we have measured apuv-, 
a uvOi an d ^break of the models in the same way we have measured them in our data. Figure 8 shows 
examples of how we measure the model spectra. The Lyman break is prominent in some models, especially 
when L/LEdd is small. However, due to relativistic boosting and aberration, small cosi (edge-on) tends to 
smear out the edges and also shift them to shorter wavelengths. Therefore, to measure chfuv for models 
with different cosz, we fit a power-law to a different smooth region of ^200A blueward of the Lyman break. 
We note that this is not exactly the same as in measuring our data, but this characterizes the UV bump of 
the models very well except for those with very strong Lyman edges (Fig. 8), and the estimate is consistent 
for all models with and without a strong Lyman edge. The NUV-optical spectral index is measured the 
same way as in our real data by fitting a power-law to two continuum windows around 1350A and 5630A 
(continuum windows b and d in Table 3). Xbreak is calculated from ccfuv an d auvo- 

Similar to our data, the distribution of auvo measured from all the models is relatively narrow with 
a median value of —1.74. We further select only models with similar black hole masses {Mbh ~ 10 8 — 4 x 
10 9 M Q ) and Eddington ratios {L/LEdd ~ 0.03 — 0.3) to those of our sample and compare them with our 
data (Fig. 7a). The auvo of the selected models has a median value of —1.91 (standard deviation a — 0.26), 
roughly in agreement with our data (median auvo = —1-52, a = 0.35). While there are no selected models 
that are as red as auvo > -1 (Fig- 7a), two objects, IRAS F07546+3928 and PG 1351+640, are redder than 
auvo > ~ 1 ! but both show obvious evidence of dust reddening. We note that 7 objects in our sample have 
L/LEdd > 0.3, and 5 objects have Mbh < 10 8 M Q , and these values of the parameters have not been covered 
by our current models. Examining the model trends in Figure 3 of Blaes (2004), extrapolating L/LEdd of the 
models to 0.5 would probably not change auvo by more than 0.1. However, lower black hole masses result 
in smaller (bluer) auvo m the models, and such values are not seen in the data. On the other hand, we 
have assumed a near maximal black hole spin in all the models used here. Lower black hole spins generally 
increase auvo m the models, and would improve agreement with the data. It would therefore be worth 
exploring such low spin models in the future. 

apuv from the models still has a broad distribution, also in agreement with our data (Fig. 7b). Since 
the distribution of auvo is narrow, the spectral breaks at the UV bump {Xbreak and Aa = apuv ~ auvo) 
are mainly defined by the change of FUV spectral index apuv in both our data and the models. aFuv is 
sensitive to Mbh, L/LEdd, and cosi, but not to the viscosity parameter a V i SCO , as shown in Figure 8. 

We further compare our data with the models by showing the changes of apuv, and Xbreak with 
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Fig. 8. — Examples of power-law fitting to the disk model spectra. Solid lines are the model spectra, and 
dotted lines are the fitted power-laws. Unless marked in each panel, the default model parameters are: 
Mbh = 1 x 10 9 M©, L/LEdd — 0.145, a V i SCO = 0.1, and cosi = 0.8 (i = 37°, close to face-on). Note that 
FUV slopes obviously change with Mbh, L/LEdd and cosz (but not a V i SCO ), but UV-optical slopes do not 
change much. Also note the strong Lyman break in the spectrum with low L/LEdd (0.0022). 
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Mbh and L/L<Edd (Fig. 9 and 10). We chose the models for extreme face-on (cosi=0.99) and extreme 
edge-on (cosi=0.01) cases with a V i SCO =0.1. 

In the face-on models, there is a clear correlation between ctFuv (hence Aq) and Mbh for L/LEdd ~ 
0.01 — 0.3 (Fig. 9, top-left). However, no evidence of such a trend can be seen in our data, even when we 
group the objects into subsamples of narrower L/LEdd ranges. Moreover, the models do not cover the otFuv- 
Mqu space with enough overlap of the data except for extremely small L/LEdd, but these small L/LEdd 
are not seen for for our AGNs. Some AGNs seems to follow the model prediction in one plot (2-dimensional 
space) , but they do not match the same model prediction in other plots (other dimensions) . Although we do 
not have information on the inclinations of objects in our sample, large disk inclinations only increase the 
discrepancy (Fig. 9, top-right). We will discuss this inconsistency more in §6. 

Xbreak does not seem to correlate with Mbh (Fig. 9, bottom-right). For large L/LEdd in face-on cases, 
Xbreak is close to 1000A, but for smaller L/LEdd, Xbreak goes to longer wavelengths. This is because olfuv 
is underestimated for models with a strong Lyman edge, as can be seen in Figure 8 (bottom-right). The 
distribution of model Xbreak also shows this clearly (Fig. 11). While there is a peak around 1000A in the 
histogram, there is also a second bump between 1200-1400A which accounts for the models with a strong 
Lyman edge. This bump disappears in the case of the selected models, where models with small L/LEdd 
(strong Lyman edge) are excluded. In the edge-on cases (Fig. 9, bottom- right), large inclination causes 
strong relativistic effects which smear out the Lyman edge, bringing the measured Xbreak back to around 
1000A. The extremely large values of Xbreak for L/LEdd=0-289 simply indicate that there is not a clear 
spectral break, because oleuv ~ o^uvo- 

We have also plotted oifuv, Aa, and Xbreak against L/LEdd for different Mbh (Fig. 10). No clear 
correlation is seen. If we could extend L/LEdd of the models to above 0.3, the models seem to cover a 
region that overlaps with our data points (face-on), but the Mbh required for the models needs to extend 
above 10 9 M Q , much higher than those calculated for most of our AGNs. The apparent correlation seen in 
Xbreak- L/LEdd is largely biased by models with small L/LEdd, for which Xbreak is overestimated due to the 
strong Lyman edge. 

Large inclination angle (small cosz, edge-on) has several effects on the models: (1) it decreases (flattens) 
oifuv, weakening the correlation between oifuv and Mbh seen in face-on models (Fig. 9, top-left); (2) it 
results in smaller difference between oifuv and osjjvo and tends to smear out the UV bump (oifuv ~ Quvo), 
resulting in unrealistically large (or small) Xbreak (Fig. 9, bottom- right); (3) it forces the models to a narrow 
region in ctpuv-L/ LEdd and Aa-L/Lsdd space regardless of Mbh (Fig. 10). 

6. DISCUSSION 

Our data set of quasi-simultaneous FUV-to-optical spectrophotometry is the first of its kind. These 
spectra, with their FUV coverage, are extremely useful for studies of AGN SEDs, especially the UV bump 
and spectral break associated with the Lyman limit. Most of our AGNs show a spectral break around 
HOOA (with large uncertainty for some objects), similar to what is seen in Zheng et al. (1997), and is in 
agreement with the non-LTE disk models. The distribution of UV-optical spectral indices redward of the 
break, and far-UV indices shortward of the break, are also in rough agreement with the models. However, 
we do not see a correlation between the far-UV spectral index and the black hole mass, as predicted by 
the face-on models. Moreover, our AGNs occupy a region in osfuv-Mbh space that is not covered by the 
thin-disk models with L/LEdd in the range of 0.01-0.3 covered by the models. These findings imply that 



cosi = 0.99 (face-on) 



cosi = 0.01 (edge — on) 




Fig. 9. — Properties of the spectral break vs. black hole mass for both data and model predictions (face-on 
and edge-on). The different open symbols are for the individual objects in different L/LEdd ranges. The 
solid lines are the model predictions. L/LEdd for each line is marked with a number close to the lines. The 
irregular patterns (intersections) of the models in the lower panels are due to the uncertainty of calculated 
Abreafc, which is mainly from the way we measure cxfuv m models with strong Lyman edge (see text). 
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Fig. 10. — Properties of the spectral break vs. Eddington ratio for both data and model predictions (face-on 
and edge-on). The different open symbols are for the individual objects in different Mbh ranges. The solid 
lines are the model predictions. Mbh for each line is marked with a number close to the lines in units of 
M 9 (M 9 = 1O 9 M ). 
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Selected Models. 




1000 2000 



Fig. 11. — Distribution of Xbreak of the sample (solid line), all models (dotted line), and the selected models 
(dashed line, 10 8 Af© <M M <4x 10 9 M Q , 0.03 < L/L Edd < 0.3). Note the peak for the models is around 
1000A. The second bump between 1200-1400A for all models is due to an overestimate for models with a 
strong Lyman edge. 



some fundamental assumptions in the models and/or in our understanding of AGN phenomena may need 
rethinking. We discuss a few relevant issues below. 

(1) We do not have information on disk inclination for our sample. Model predictions show that large 
inclinations do change ckfuv, but comparing the edge-on and face-on cases, a large inclination increases 
the disagreement between our data and the model predictions in the parameter space (a fuv j Xbreak j Mbh , 

L/LEdd)- 

(2) Reddening can increase the observed apuv- We note that our data will match the models better 
in olfuv-Mbh space if all the a fuv are shifted down by ~1. This will also allow us to have non- face-on 
inclinations for our sample, and it is unlikely that we systematically under-corrected the Galactic redden- 
ing. Scott et al. (2004) looked for possible systematic effects that could be caused by incorrect reddening 
corrections in the FUSE observations of low-redshift AGN and found none. However, any intrinsic redden- 
ing can play an important role. Two objects in our sample, IRAS F07546+3928 and PG 1351+640 show 
evidence of intrinsic reddening. Both of them show intrinsic absorption features, and unlike other objects, 
their FUV to blue optical spectra significantly deviate from a power-law (see Fig. 1). All the above evidence 
suggests a strong (intrinsic) reddening effect. If we assume the intrinsic reddening has the same nature as 
the Galactic reddening and follows the same extinction law (Cardelli et al. 1989, CCM), roughly a correction 
of E(B-V) = 0.03 is needed to bring apuv down by about 1; for the extinction curve of the Small Magellanic 
Cloud (SMC, Prevot et al. 1984), this requires E(B-V) = 0.04. These are modest amounts, but they have 
significant effects at short wavelengths. 

To show the effect of reddening on the spectral break, we performed some simple simulations (Fig. 12). 
We found that for the CCM reddening curve (Fig. 12a), adding reddening of E(B-V) = 0.04 to a power- 
law spectrum can produce a UV turnover that resembles the spectral break we see in AGN spectra. Also 
dcrcddcning a spectrum with a spectral break can virtually eliminate the break. Compared to the intrinsic 
reddening seen in some nearby AGNs (Maiolino et al. 2001), the reddening values required here are very low 
and would not be surprising to find in AGNs. On the other hand, if we use the flatter SMC extinction curve, 
a larger E(B-V) (> 0.04) is needed to produce noticeable results (Fig. 12b). However, it cannot produce a 
clean UV turnover, or remove a spectral break without introducing a large-scale curvature to the spectrum 
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Fig. 12. Simulations of reddening a power-law continuum {auvo = —1-52, the median value of our sample), 
and dereddening a broken power-law continuum {auvo = —1-52, aFuv=0)- ( a ) reddening a power-law (top, 
thick line) and dereddening a broken-power (bottom, thick line) with the Galactic reddening curve (Cardclli 
ct al. 1989). Note the spectral break at 1000A disappears (bottom) when E(B-V)=0.04. (b) same as (a) 
with the SMC reddening curve (Prevot et al. 1984). Note the curvature between 800-3000 A (top) for 
E(B-V)=0.08. 

in the NUV to optical region. IRAS F07546+3928 and PG 1351+640 in our sample seem to show this 
curvature (Fig 1). 

We note that the extinction curves we use above for the FUV region are simple extrapolations from 
CCM and SMC extinction curves. Hutchings & Giasson (2001) derived FUV extinction curves for stars in 
the Galaxy, Large Magellanic Cloud (LMC), and SMC using FUSE data, and found that they appear to 
extend the extinction curves from longer wavelengths in a straightforward way. Sasseen et al. (2002) also 
found an FUV extinction curve in the Galactic diffuse interstellar medium consistent with an extrapolation 
of the CCM curve (Ry = 3.1), but it is likely that the Galactic extinction curve is not applicable to AGNs, 
and we see significant difference between the CCM and SMC curves in the simulations. In fact, Maiolino 
et al. (2001) found that the ratio of E(B-V) to hydrogen column density Nhi in AGNs is lower than the 
Galactic value by a factor of <~ 3 — 100, and suggested that the dust in the circumnuclear region of AGNs has 
different properties than in the Galactic diffuse interstellar medium. In a study of red and reddened quasars 
in SDSS, Richards et al. (2003) found that an SMC-like reddening law with E(B-V) between 0.135 and 0.07 
can redden their normal color composites to a dust-reddened composite spectrum. Normal quasars in the 
SDSS exhibit little to no intrinsic reddening. Hopkins et al. (2004) find that 81% of the SDSS quasars have 
E(B-V) < 0.02, and those that are reddened follow an SMC-like extinction law. With a simple assumption 
that all AGNs have the same continuum slope, Gaskell et al. (2004) derived an extinction curve for AGNs, 
and claimed that it is much flatter than the Galactic CCM extinction curve. Whether this is true or not, if 
the reddening curve in the AGNs is at least as flat as the SMC curve, the reddening is not able to produce 
the spectral break seen in our AGNs without leaving a clear signature at longer wavelengths. In addition, as 
seen in Figure 12b for the SMC law, a relatively small reddening (E(B-V) < 0.08) can significantly suppress 
the UV continuum, as is only seen in a few of our objects. This implies that either the intrinsic reddening 
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for most objects in our sample is very low, or the intrinsic reddening curve for AGNs is very different (e.g., 
flatter) from the SMC curve. 

It should be kept in mind that these AGNs were targeted for observation by FUSE because they were 
known to be bright in the UV. This selection would bias our sample to be among the AGNs with the least 
intrinsic reddening. If our results indeed arise from this effect, they may be stronger and more common 
among the general AGN population than in our sample. 

Until reliable extinction curves and quantitative intrinsic reddening for individual AGNs are available, 
reddening effects and the AGN intrinsic continuum slope cannot be completely decoupled. As discussed 
above, intrinsic dust reddening could significantly modulate the AGN UV bump, but it cannot produce a 
clear spectral break if the extinction curve is flatter than the CCM law. The observed spectral break is 
indeed intrinsic to AGNs. 

(3) Comptonization can also alter the UV and soft X-ray spectral indices as well as smear out the Lyman 
edge (e.g., Czerny & Zbyszewska 1991; Hubeny et al. 2001). Zheng et al. (1997) were able to fit their HST 
composite spectrum with a disk model plus Comptonization (Czerny & Zbyszewska 1991); Kriss et al. (1999) 
did the same for 3C273. With the same SED of 3C273, Blaes et al. (2001) found that while Comptonization 
has no effect on the optical and mid-UV spectrum, where they got a reasonable fit with the non-LTE disk 
model by Hubeny et al. (2000), including Comptonization is necessary to smear out the Lyman edge feature 
and to extend the disk spectrum to the soft X-ray band. It may also be true that Comptonization is 
important only in some objects, but the models we studied do not included the Comptonization in creating 
the integrated spectra. Adding Comptonization to disk models will introduce a scattering medium and hence 
more parameters. We also note that the effect of Comptonization on the model spectra is very similar to 
the relativistic smearing effect which is strong in large inclination disks. 

(4) We use geometrically thin disk models to compare with our data. The L/LEdd for many of our 
objects exceed the thin disk model limit of 0.3, and there is no physical reason why L/LEdd cannot exceed 
this limit. Increasing L/LEdd would be expected to transform a thin disk into a slim disk (Szuszkicwicz, 
Malkan, & Abramowicz 1996), for which advection may become important (Blaes et al. 2001). Not only 
does a slim disk model produce L/Lsdd that is consistent with values estimated for many quasars, it could 
also improve the fit to observed spectra as suggested by Blaes et al. (2001). Models with nonzero magnetic 
torques across the innermost stable circular orbit (Agol & Krolik 2000) may also improve the fit to observed 
spectra. Detailed model spectra are needed to compare with our data. 

(5) We used models with a fixed value of a V i SCO and maximum black hole spin to compare to our UV 
data. While the FUV slope is not sensitive to a V i SCO , if our objects indeed have very different black hole 
spins, the predicted correlation between qfuv and Mbh for a single value of spin may not be seen in the 
data. In our future work we will construct specific models for the masses and luminosities measured for our 
objects, and try to find for each object the best-fit inclination and black hole spin. Inclinations can also be 
roughly constrained using the radio properties of radio- loud AGNs (e.g., Orr & Browne 1982). With fewer 
free parameters, the models can be better constrained and tested by observational data. 

7. SUMMARY 

1. We construct SEDs of 17 AGNs with quasi-simultaneous spectra covering 900-9000A (rest frame). The 
SEDs are available in digital format at http://physics.uwyo.edu/agn/. 
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2. The distribution of ajjvo is narrow, and in rough agreement with non-LTE thin-disk models. The 
distribution of cxfuv of our sample is also in rough agreement with that of the models. 

3. We see a spectral break in the UV for most of our objects, and the break is around llOOA. Although 
this result is formally associated with large uncertainty for some objects, the FUV spectral region 
is below the extrapolation of the NUV-optical slope, indicating a spectral break around llOOA, in 
agreement with previous studies of HST composite spectra. 

4. Intrinsic dust reddening can significantly modulate the AGN continua, but the spectral break is intrinsic 
to the AGNs, and is not caused by possible reddening if the dust extinction curve in AGNs is flatter 
than the Galactic reddening curve. 

5. We do not find the correlation between apuv an d Mbh expected by the thin accretion (face-on) disk 
model, possibly due to the small sample size. Scatter introduced by other varying disk parameters 
that are not included in our models, such as inclination and the black hole spin, could also weaken the 
expected correlation. 

6. Thin-disk models do not match the observed spectra in the space of otpuv an d Mbh within the thin- 
disk model limit (L/LBdd=0.3). This discrepancy may be attributable to the effects of Comptonization 
and other factors the models have not included. 

We thank Beverley Wills for providing the optical spectra of PG1100+770. This work is based on data 
obtained for the Guaranteed Time Team by the NASA-CNES-CSA FUSE mission operated by the Johns 
Hopkins University. This work has been supported by NASA through grants from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under 
NASA contracts NAS5-32985 and NAS5-26555. 

Facilities: FUSE, HST(STIS), KPNO. 
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